Amplitude-mode dynamics of polariton condensates 
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We study the stability of collective amplitude excitations in non-equilibrium polariton conden- 
sates. These excitations correspond to renormalized upper polaritons and to the collective amplitude 
modes of atomic gases and superconductors. They would be present following a quantum quench or 
could be created directly by resonant excitation. We show that uniform amplitude excitations are 
unstable to the production of excitations at finite wavevectors, leading to the formation of density- 
modulated phases. The physical processes causing the instabilities can be understood by analogy 
to optical parametric oscillators and the atomic Bose supernova. 



In addition to the phase mode responsible for su- 
perconductivity and superfluidity, Bose-Einstein conden- 
sates possess collective amplitude modes. Examples are 
the amplitude modes of superconductors [T] and cold 
atomic gases |lHH]j in which the density fluctuates be- 
tween condensed and non-condensed particles, and the 
formally equivalent Higgs mode of a relativistic scalar 
condensate [9l . These modes are orthogonal to the phase 
modes in the order-parameter space, and their excitation 
is predicted to lead to non-equilibrium states in which the 
magnitude of the order parameter oscillates [2HS1 IHl IH] ■ 

The strong coupling of excitons and photons in semi- 
conductor microcavities leads to the formation of upper 
and lower polaritons. Pumping such a microcavity is 
observed to lead to the formation of Bose-Einstein con- 
densed states of polaritons, characterized by coherent 
light emission from the structure |12j . This system al- 
lows the investigation of physics analogous to that of 
cold atoms, but in a strongly coupled regime with long 
range interactions mediated by the photons. The cou- 
pling of polaritons to electromagnetic radiation outside 
the microcavity gives an advantage over cold atoms in 
that polaritons can be coherently controlled by external 
pumping, and conversely their dynamics and coherence 
can be be directly observed in the emitted light [12^ JJJ . 
This has allowed experiments to reveal the collective be- 
havior arising from excitations of the phase mode, includ- 
ing superfluidity [14 , vortex dynamics jl5H17j . and the 
Bogoliubov spectrum [TSl [TH]. However, the collective 
behavior associated with the amplitude mode has yet 
to be considered. Controlled excitation of the amplitude 
mode of an atomic Fermi gas requires a rapid switch of 
the magnetic field [IHH]) while in superconductors a co- 
existing charge-density wave is needed [I]. 

The amplitude mode of the polariton condensate may 
be identified by considering the Dicke model [501 1^ . 
which for weak coupling becomes the BCS model. This 
limit, with well-known amplitude and phase modes, is 
adiabatically connected to the strong-coupling limit re- 
alized in microcavity experiments, where the collective 



modes are upper and lower polaritons. Thus the ampli- 
tude mode should be identified with the upper polari- 
ton. Uniquely, in this system the amplitude mode can 
be directly driven by resonant excitation, allowing the 
resulting collective behavior to be studied experimen- 
tally. In this paper we predict the collective behavior 
arising when the amplitude degree-of-freedom of a po- 
lariton condensate is manipulated in this way. Whereas 
driving the phase mode induces superflows, we find that 
driving the amplitude mode completely destabilizes the 
condensate, causing the polaritons to spontaneously or- 
ganize into density-modulated phases (Fig. [2]). This oc- 
curs because the interactions transfer the excess rest mass 
and interaction energy of the non-equilibrium state into 
kinetic energy, in a way ruled out for the phase modes 
by the Landau criterion. 

We begin by considering the dynamics of a Dicke model 
of polariton condensation, because for this model exact 
solutions for the collective dynamics are available [6]- 
|S1 [22] • We show that these solutions, which describe 
a uniform oscillating condensate, are unstable once exci- 
tations at finite wavevector are considered, and system- 
atically identify the instabilities. Our analysis reveals 
two types (Fig. [T]): a wave- mixing instability between 
the lower and upper polariton, and a modulational in- 
stability caused by an attractive interaction between up- 
per polaritons. We propose a Ginzburg-Landau theory 
which captures these instabilities, and provides a realis- 
tic model of a microcavity. We use this theory to predict 
the true steady-states under continuous excitation. 

We can systematically establish the amplitude-mode 
dynamics of a polariton condensate by considering the 
generalized Dicke model [501 [H] ■ 

i,k 

This model describes N localized exciton states with po- 
sitions r^, in the limit where exciton-exciton interactions 
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exclude double occupancy. Thus each state may be oc- 
cupied (erf = 1) or unoccupied (erf = — 1). The excitons 
are coupled to the two-dimensional microcavity photons 
with in-plane wavevectors fc, annihilation operators "0^, 
and dispersion relation cok ~ a;o+|A:p/2mph (^-=1), where 
mph ~ lO^^mc- At this stage we are neglecting impor- 
tant effects such as the dispersion and inhomogeneous 
broadening of the exciton, the polarization degrees-of- 
freedom [321 Hll : and the finite lifetime of the cavity pho- 
tons. We consider only the limit of a single exciton en- 
ergy E, and a single coupling strength parametrized by 
the Rabi splitting flu. To treat the collective dynamics of 
([T]) it is convenient to work with the expectation values 

= {A)/VN, Pk = {Pk)/VN = ^E,('^r>e-^'"'% 
and Dk = {bk)/^ = F (-^f) e"'''"'- and Pk are 
the macroscopic components of the electric field and po- 
larization at wavevector fc, while measures the exciton 
occupation. 

With only a single photon mode the Heisenberg equa- 
tions for the model ([T]) are integrable, and closed-form 
solutions for the dynamics of the collective variables are 
available [BHSl [22] • These solutions have been studied 
in the context of quench experiments on atomic Fermi 
gases, in which the pairing interaction would be rapidly 
switched, giving an impulsive excitation of the collective 
amplitude mode (analogous experiments are proposed in 
light-matter systems 10, llj). This is predicted to lead 
to a spatially uniform condensate in which the order pa- 
rameter oscillates in time, of which an example is jT^ : 

V'o - 0(i) = 0+dn((/.+t,At), (2) 

Po = p{t) = 2[-wo<^(t) + i^m/rtR, (3) 

where dn is a Jacobi elliptic function, and the zero of 
energy is such that V'o is real. The parameters cj)^ and k, 
determine the period T and magnitude of the oscillations, 
and are known functions of the model parameters and 
the initial conditions. Note that the oscillation frequency 

= 27r/r generally differs from due to interactions. 

Since a microcavity supports a continuum of in-plane 
modes we must consider the behavior of an oscillating 
condensate beyond this mean-field approximation. To 
do this we linearize the equations of motion about the 
mean-field solution ijjQ = 4){t), Pq = p{t) and Z?o = ^(^)- 
The fluctuating parts of the collective variables obey 

iS-^k = i^k^k + flRSPk/2, 
i5Pk = ESPk - r!fl(A<5Vfc + ^SDk)/2, 
idbk = {p*S^Jk + ^SPlk - 4>*5Pk - pSr^k) , (4) 

where the coefficients are time dependent because the 
condensate oscillates with angular frequency VL, which 
is the gap energy to the occupied amplitude mode. In 
deriving Q we assume that the wavevectors are much 
smaller than the inverse of the spacing of the exciton 
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FIG. 1: (Color online). Quasi-energy spectrum of an oscil- 
lating polariton condensate. Lines (left axis) correspond to 
the real part of the quasi-energy A, and blue crosses to the 
non-zero imaginary parts (right axis). Modes with non-zero 
£sA are unstable. The line colors indicate modes deriving 
from the phase mode/lower polariton (red solid), or the am- 
plitude mode/upper polariton (black dashed). Inset: spec- 
trum of a weakly-interacting Bose-Einstein condensate with 
attractive interactions, qualitatively reproducing the small- 
wavevector behavior. The parameters are Wo = —E — OA^lji, 
A(0) « 0.71, (j){0) « 0.44, implying a renormalized oscillation 
frequency fl ~ 0.480^. 

states, SO that motional narrowing is effective and the 
wavevector of the fluctuations is well defined [25j [26] . 

Eqs. ([4]) have the form ^t* = A(t)'^ where is a vector 
of the fluctuating fields and A{t) a time-periodic matrix. 
Thus Floquet's theorem applies and the solutions are of 
the form ^ = u{t)e'^^*, where u{t + T) = u{t) and the 
quasi-energy A is defined up to an integer multiple of 
fl. Since u is periodic, the stability of the condensate is 
determined by the imaginary part of A. 

To illustrate the fluctuation spectrum of an oscillating 
condensate, we show in Fig. [ij the result obtained for the 
exact solution (2]3). We measure energy and frequency 
in units of the Rabi splitting Jl/j, and wavevector in units 
of ko = y/2mphilB. The spectrum can be understood by 
considering the limit of small oscillations around an equi- 
librium condensate, where ([2| becomes (j) — (f>o+(j)i cos fit 
with 01 ^ 00 1 and we have taken the zero of energy 
to be the equilibrium chemical potential. In this limit 
the construction of the spectrum is analogous to that of 
the nearly free electron model |27j: with the weak time- 
periodic component playing the role of the weak periodic 
potential. To lowest order in the oscillating component, 
01 = 0, the fluctuation spectrum is that of an equilibrium 
polariton condensate [20l [2Tj . There is a mode starting 
at zero energy and wavevector, which in the low-density 
limit is the lower polariton, and in Fig. [T] has been renor- 
malized into the linearly-dispersing phase mode. There 
is also a gapped mode at zero wavevector, which in the 
low-density limit is the upper polariton [21], and in Fig.[l] 
is the collective amplitude mode appearing at the gap fre- 
quency [2H] • The effect of the oscillations is to fold the 
spectrum in frequency, and to couple together the result- 
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ing levels. Depending on the phase of the coupling this 
can result in either a level repulsion or attraction. In the 
latter case the dispersion relation is flattened, and imagi- 
nary parts appear for the quasi-energies. This signals an 
instability of the spatially uniform solution, and an ini- 
tially exponential growth of modes at finite wavevectors. 

The two strongest instabilities in Fig. [l] occur near 
|fc| = |fci| ~ O.Bfcp and \k\ « 0. The first occurs where 
the positive-energy branch derived from the lower polari- 
ton crosses with a replica of the corresponding negative- 
energy branch. In the low-density limit this occurs only 
for E ^ ujQ > 0. The result is a wave- mixing instability 
in which an upper and lower polariton from the oscil- 
lating condensate at fc = scatter to a pair of lower 
polaritons at ±ki. This is analogous to the instability 
that drives the microcavity optical parametric oscillator 
(OPO) [29 , with the two k = states forming the pump, 
and the states at ±fci the signal and idler. This scattering 
process has previously been considered as a loss mecha- 
nism for incoherent polaritons [30, . The creation of phase 
modes from amplitude oscillations has been considered in 
the Bose-Hubbard model [3] and a similar instability has 
recently been found in the BCS model ,31j. 

The second instability in Fig. [T] corresponds to the flat 
dispersion relation in the upper polariton branch near 
k = 0. It occurs because the saturation of the light- 
matter coupling [32] reduces the Rabi splitting with in- 
creasing excitation. Thus there is an attractive inter- 
action between upper polaritons, and states containing 
more than a single such excitation are unstable. The re- 
sulting form of unstable spectrum is that obtained from 
the Bogoliubov analysis for a condensate with weak at- 
tractive interactions (inset). 

The microscopic theory is capable of describing the 
initial instability, but becomes unwieldy in the nonlinear 
regime as the unstable modes evolve. To consider the 
long-term evolution in a realistic microcavity we study 
the Ginzburg-Landau theory 

-i"f^p + ^ + F, (5) 

Here P{x,t){tfj[x,t]) represents the macroscopically- 
occupied exciton (photon) field, which is linearly cou- 
pled to the photon (exciton) field to generate polaritons. 
The nonlinear coupling accounts for the saturation ef- 
fect. This form can be obtained from ([I]) by represent- 
ing the exciton operators using the Holstein-PrimakofF 
transformation, and follows directly from the microcav- 
ity exciton-photon Hamiltonian |32j treated in mean-field 
theory. We neglect the Coulomb interactions between the 
excitons because they act within each quantum well, and 
hence are weaker than the saturation nonlinearity when 
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FIG. 2: (Color online). Photon density from numerical sim- 
ulations of Eq. ([5|, with resonant pumping of both upper 
and lower polaritons (left panel), and the upper polariton 
alone (right panel). The left panel shows the formation of a 
short-wavelength density modulation due to the wave-mixing 
instability, and the right the result of the attractive-gas in- 
stability. The pump field profile is a Gaussian at a; = with 
CT = 10 /im; D,R = 20 meV; E — luq = A meV. 

the number of wells is large. This is consistent with the 
observed redshifting of the upper polariton with density 
in some microcavities [33] . We have introduced a damp- 
ing constant 7 — 0.25 ps~^ to describe decay of the cav- 
ity photons, a term F to describe a resonant pump laser, 
and a noise source ^ to model spontaneous emission noise. 
We focus on resonant excitation with circularly-polarized 
light, and hence include only a single polarization of ex- 
citon and photon. For numerical work we set A = 1, 
absorbing the nonlinear coupling strength into the defi- 
nition of density. 

Fig.[2]shows the photon densities calculated from ^ in 
one space dimension, with resonant excitation of both the 
upper and lower polariton (left panel), and the upper po- 
lariton alone (right panel). For these parameters pump- 
ing both modes initially leads to a uniform induced con- 
densate within the pump spot, in which as in ^ the field 
amplitude oscillates due to intermode beating. However, 
at later times this uniform state breaks down, and we see 
the formation of a condensate with a short-wavelength 
density modulation. We do not see strong signs of ther- 
malization and chaos developing from the instability, as 
has been suggested for atomic systems [31 |3J , presum- 
ably due to the presence here of dissipation. 

In Fig. [2] we also show the behavior when only the 
upper polariton is pumped. In this case the wave-mixing 
instability cannot occur, and the dominant instability is 
due to the attractive interactions. In cold atomic gases, 
producing a condensate with attractive interactions leads 
to the Bose supernova [31] where the condensate explodes 
due to the excess interaction energy of the uniform state. 
In the microcavity, we instead predict the formation of 
a large wavelength density modulation. This is because 
the polaritons can organize in such a way that the excess 
energy injected by the pump is dissipated, leading to a 
relatively stable steady-state. 

While in the undamped model (|l|) the uniform state is 
always unstable, the dissipation inTsj) implies a threshold 
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density below which uniform states are stable. This den- 
sity no is where the gain due to the interactions exceeds 
the losses from the modes, and therefore for an inter- 
action strength g ~ ^rX is guQ ~ 7, up to numerical 
factors which are typically of order one. This is essen- 
tially the threshold criterion for the OPO [53], so that the 
threshold density should be within reach experimentally. 

The models we have considered can be extended to 
include more realistic details. In particular, we have 
neglected the electron-hole continuum, which can pro- 
duce resonant damping of the upper polariton, similar 
to the amplitude-mode damping predicted in atomic sys- 
tems [28j [35l [36] . However, the two-dimensional exciton 
binding energy in GaAs (CdTe) is 20 (40) meV, so that 
the upper polariton in Fig. [2] is below the continuum, and 
not too strongly damped by this mechanism. 

An interesting extension of our work would be to al- 
low both polarizations of polaritons and excitons in the 
Ginzburg-Landau theory. In this case there will be addi- 
tional amplitude modes connected to fluctuations in the 
degree of polarization of the condensate. It would be 
interesting to determine whether the polarization oscilla- 
tions previously seen in OPO simulations |23j correspond 
to these amplitude modes, and to investigate the possibil- 
ity of nonlinear decay processes similar to those discussed 
here. It would also be useful to investigate the possibil- 
ity of attractive interactions between amplitude modes in 
the BCS model, and hence establish the extent to which 
the instabilities identified here occur in other systems. 

In summary, we have used the Dicke model to show 
that a uniform condensate in which the amplitude mode 
is excited is unstable due to (a)an attractive interaction 
between amplitude modes and (b) scattering between am- 
plitude and phase modes. We have used a Ginzburg- 
Landau theory to show that these instabilities lead to 
the formation of spatially inhomogeneous condensates. 
In a microcavity the amplitude mode corresponds to the 
upper polariton, and therefore these instabilities can be 
induced by resonant excitation, leading to features in the 
real-space density and to bright emission at an angle from 
the cavity. Our work shows that there is a rich collective 
behavior associated with condensate amplitude modes, 
and that microcavities provide a unique opportunity to 
explore this physics experimentally. 
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